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Abstract 

We construct a local interacting quantum dimer model on the square lattice, whose zero- 
temperature phase diagram is characterized by a line of critical points separating two or- 
dered phases of the valence bond crystal type. On one side, the line of critical points ter- 
minates in a quantum transition inherited from a Kosterlitz-Thouless transition in an asso- 
ciated classical model. We also discuss the effect of a longer-range dimer interaction that 
can be used to suppress the line of critical points by gradually shrinking it to a single point. 
Finally, we propose a way to generalize the quantum Hamiltonian to a dilute dimer model 
in presence of monomers and we qualitatively discuss the phase diagram. 
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Introduction 



Dimer models are of interest to a variety of scientific disciplines from chemistry 
to mathematics and physics. In chemistry, dimers are used, for example, to model 
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molecules deposited on crystalline surfaces and to study their thermodynamic prop- 
erties [T]. In mathematics, dimers are often used to construct combinatorial and 
folding problems such as the domino tiling of a two-dimensional (2D) plane [2]. 
In physics, dimer models have been elevated from problems in classical statistical 
physics itl , H, 5 B, B, [ sj I9 , 10], to problems in quantum statistical physics ITT 
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J], with the advent of high-T^. superconductivity. In particu- 
lar, quantum dimer models can provide examples of strongly correlated quantum 
systems for which the zero temperature phase diagram is characterized by exotic 
quantum phase transitions that fall out of the classification of phase transitions pro- 
posed by Landau [|40l. kill. Furthermore, the finite-temperature phase diagram of 
quantum dimer models might give some insight into the phenomenological obser- 
vation that scaling laws extend to surprisingly high temperatures in some strongly 



correlated systems [4^ |43|] 



In this paper, we show how dimer models can be used as a laboratory to construct 
quantum Hamiltonians displaying phase transitions that cannot be understood in 
terms of a local order parameter, i.e., phase transitions that cannot be encoded by 
an effective theory of the Landau-Ginzburg type, a topic of renewed interest in 
condensed matter physics [|40l.l4lll. Perhaps the most famous counter example to a 
phase transition described with a Landau-Ginzburg action for a local order parame- 
ter is the Kosterlitz-Thouless (KT) transition. The KT transition is a weak essential 
singularity of the free energy for a phase-like order parameter with support in 2D 
Euclidean space. It is interpreted as the unbinding of topological defects (vortices) 
in the order parameter. The main result of this paper is the construction of a quan- 
tum dimer model with local interactions on the square lattice [Eqs. (13,14)] that 
undergoes a quantum phase transition of the KT type when measured by the spatial 
decay of equal-time correlation functions. 

It is well known that quantum phase transitions can be of the KT type in ID sys- 
tems with dynamical exponent z = 1 relating the scaling in space to the scaling 
in time. For example, a ID Luttinger liquid can be unstable to a charge-ordered 
density wave through a KT transition. This is so because the quantum field theory 
describing the quantum phase transition for interacting fermions is related through 
bosonization to a scalar field theory, the Sine-Gordon model. Analytical continu- 
ation of time to imaginary time can be used to turn (Minkowski) space-time into 
2D Euclidean space while the Poincare symmetry group becomes symmetry under 
translations and rotations. Evidently, if a quantum phase transition can be described 
by a local quantum field theory in D + 1 space-time that turns into a local classical 
action undergoing a classical phase transition inD+l Euclidean space upon analyt- 
ical continuation of time to imaginary time ll44n . then this quantum phase transition 
cannot be associated to a KT transition when z = 1 and D > 2. Viewed against 
this no-go theorem, it is remarkable that some equal-time correlation functions of 
a 2D quantum dimer model with local interactions and defined on the square lat- 
tice [Eqs. (13,14)] share the hallmarks of a KT transition. 
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Since the quantum Hamiltonian of our 2D lattice model, defined in Eqs. (13) 
and (14), has only local interactions, it is safe to argue that unequal-time corre- 
lations should show algebraic behavior if the equal-time correlations do so. The 
reason is that a local Hamiltonian with algebraic spatial correlations should be gap- 
less. A rigorous proof of this statement was given by Hastings in Ref. ll45ll . The 
converse statement is known not to be true, as shown in Ref. ll46ll . A local quan- 
tum Hamiltonian may be gapless but have only short-ranged spatial correlations 
between local operators. Therefore, while in this paper we concentrate solely on 
equal-time correlation functions, the KT-like transition identified through the alge- 
braic spatial correlations should be manifest in unequal-time correlation functions 
related to the spatial ones through a dynamical exponent z. Moreover, a dynamical 
exponent ^ = 2 is suggested by our mapping of the quantum system onto a classical 
model with local stochastic dynamics [^16, 28, 33ll . 



What is the continuum imaginary-time field theory that captures the low-energy 
physics of the local quantum lattice Hamiltonian, Eqs. (13,14), when fine-tuned 
to its line of critical points? We expect it to be that of a classical Lifshitz point 
problem in a uniform magnetic field introduced by Grinstein in Ref. [47] . Here, a 
short-ranged anisotropic coupling between 2D layers, each of which are described 
by a local classical Lagrangian, is interpreted as the coupling between imaginary- 
time slices of the quantum problem 11291 14811 . Note, however, that it is not always 
the case that a local classical Lagrangian in D + 1 Euclidean space corresponds 
to a local quantum Hamiltonian in D spatial dimensions; a counter-example was 
shown in Ref. [|49ll . Our 2D lattice realization of the quantum KT-transition starts 
directly from the local quantum Hamiltonian, and avoids any discussion of the cor- 
responding classical imaginary-time Lagrangian. 



The quantum dimer model on the square lattice discussed in the present paper, 
Eqs. (13,14), is represented by a symmetric and positive matrix that obeys the so- 
called Stochastic Matrix Form (SMF) decomposition 11281 1291. 13311 . The advantages 
of an SMF decomposition of a quantum Hamiltonian are three-fold. First, at least 
one ground state (GS) can be obtained exactly in terms of the parameters entering 
the SMF quantum Hamiltonian, see Eq. (16) lll2lll6l l28l. 1291 13311. Second, one can 
construct a classical configuration space that is in one-to-one correspondence with 
the orthonormal basis in which the SMF quantum Hamiltonian is represented. On 
this configuration space, a classical partition function can be uniquely defined from 
the GS wavefunction, see Eq. (19) lll2l ll6l. l28U29l l33ll. such that zero-temperature 
and equal-time correlation functions of quantum operators diagonal in the SMF 
basis are equivalent to equilibrium thermal averages of corresponding quantities 
in the classical system. Note that the possibility to use classical numerical tech- 
niques, such as Monte Carlo simulations and transfer matrix calculations, in the 
same number of dimensions gives access to much larger system sizes than quantum 
techniques, such as quantum Monte Carlo or exact diagonalization routines, do. 
Third, the parameters entering the SMF quantum Hamiltonian allow us to define 
in a unique way the approach to equilibrium of the associated classical system, see 
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Refs. [I16I. I28L 13311. i.e., in a way that zero-temperature, imaginary-time correlation 
functions of operators diagonal in the SMF basis can be obtained from real-time 
correlation functions in the stochastic classical system. In particular, if the partition 
function of the associated classical system undergoes a KT transition upon vary- 
ing the quantum parameters entering the SMF quantum Hamiltonian, so does the 
equal-time GS expectation value of operators diagonal in the SMF basis. We can 
now understand how it is possible to circumvent the no-go theorem. The no-go the- 
orem assumes that a classical phase transition faithfully represents all correlation 
functions in a quantum phase transition with z = 1. In this paper, only equal-time 
GS expectation values of operators diagonal in the SMF basis are faithfully repre- 
sented by correlation functions at the KT critical point since the value of 2; 7^ 1 is 
not known rigorously. Similar results have been announced by Papanikolaou et al. 
inRef LSQl. 

The paper is organized as follows. We will show in Sec. 2 that the GS (16) 

of quantum Hamiltonian (13) defines the classical partition function (19) for inter- 
acting dimers on the square lattice 



Z{T/u) := E e-^c 



c 



Here iSq is the set of all possible classical dimer configurations on the square lattice, 
u and T are two real parameters, and A^"^ is the number of plaquettes having two 
parallel dimers in configuration C. The classical partition function (19) was studied 
numerically by Alet et al. in Ref. [10] for one sign of the interaction between the 
dimers. In Sec. 3 we extend the numerical study by Alet et al. to the other sign of 
the interaction between the dimers. The temperature of the classical partition func- 
tion (19) plays the role of a quantum coupling in the quantum Hamiltonian (13). 
Ground-state equal-time expectation values of operators diagonal in the dimer ba- 
sis are thus inherited from the thermodynamics of the classical partition function. 
In the high-temperature regime, the associated classical system exhibits a line of 
critical points. As the temperature is lowered, the classical system undergoes either 
a first-order or KT transition depending on the sign of the interactions between the 
dimers. To study the robustness of the line of critical points in the zero-temperature 
phase diagram of quantum Hamiltonian (13), we extend the range of the dimer 
interactions in Sec. 4. We show that, for one sign of the longer-range dimer inter- 
action, the line of critical points shrinks continuously upon increasing the strength 
of this longer-range interaction. Section 5 is devoted to another kind of perturba- 
tion to our quantum Hamiltonian (13): the presence of defects represented by sites 
not occupied by a dimer (monomers). We define the more general SMF quantum 
Hamiltonian (55) that accounts for the existence of monomers in a dilute dimer 
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model. Once again the GS (56) can be computed exactly 

ce5 

where, in addition to the quantities defined above, S is the set of all classical dilute 
dimer configurations on the square lattice (i.e., where each site belongs to at most 
one dimer), /i is the chemical potential for monomers, and M^. is the total number 
of monomers in configuration C. We can then establish a correspondence between 
our quantum SMF Hamiltonian and a classical dilute dimer model described by the 
partition function (57) 

Z{T/u,fi/T) := 2^ exp — = exp . 

ce5 \ ^ / Ces \ / 

This is a promising result as equal-time correlation functions for monomers can 
thereby be studied, for example, using classical Monte Carlo algorithms [id]. Of 
course, the computation of unequal-time correlation functions for monomers still 
requires the use of quantum Monte Carlo simulations or exact diagonalization tech- 
niques and is therefore limited to smaller system sizes. 



2 A square lattice interacting quantum dimer model with solvable ground 
state 



We begin our construction of a quantum square lattice dimer model (SLDM) whose 
ground state (GS) is exactly solvable by revisiting the quantum dimer model that 
was introduced by Rokhsar and Kivelson (RK) in Ref. [12]. The Hilbert space is 
given by the span of the orthonormal basis states \C) labeled by all the classical 
dimer configurations C of the square lattice. A classical dimer configuration on the 
square lattice is obtained by covering all the bonds connecting nearest-neighbor 
sites with dimers in such a way that each site is the end point of one and only one 
dimer. The set off all allowed dimer coverings will be denoted by Sq. 



The RK Hamiltonian acting on this Hilbert space is then commonly written as 
HnK = E\^ (\U){U\p+\n){U\r)-t(\U){U\r>+\n){U 



(1) 



where the summation is over all plaquettes p of the square lattice. The operators 
I rn ) ( nn Ip l n ) ( D denote projection operators onto the subspace of states 
associated with configurations C that contain two parallel dimers (vertical or hor- 
izontal, respectively) on plaquette p. The operators | III ) ( D |p and | Q ) ( III |p are 
the plaquette p flipping operators that maps any state \C) with two horizontal or 
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vertical dimers at plaquette p onto the state \C) obtained by rotating the two paral- 
lel dimers at plaquette p by 90 degrees, while it annihilates state \C) otherwise (i.e., 
if the configuration C does not contain two horizontal or vertical parallel dimers at 
plaquette p). Given any classical dimer configuration C, we shall call any plaquette 
occupied by two parallel dimers a flippable plaquette. 

The representation (1) of the RK Hamiltonian is explicitly local as no reference is 
made to the dimer covering away from plaquette p in the definitions of operators 
I n ) ( n |p, I O ) ( O |p and I □ ) ( O |p, I O ) ( □ |p- The nature of the GS of the RK 
Hamiltonian depends on the dimensionless ratio between the characteristic energies 
f G M and t > 0. When v/t <^ —1, the GS is expected to display columnar ordering 
of the dimers, i.e., each dimer has precisely two parallel neighboring dimers. When 
v/t ^ +1, the GS is expected to display a staggered ordering of the dimers, i.e., 
no two parallel neighboring dimers are present in the system. The so-called RK 
(critical) point 

v = t (2) 



in parameter space is special in that the GS is known exactly [|12(1 



I*) := E |C) (3) 



and it is non-degenerate within each irreducible sector of the Hilbert space under 



the action of the RK Hamiltonian. It separates two ordered phases [|36|1 . 



The notation used so far is not suitable to construct generalizations of the RK 
Hamiltonian. What we need instead is a representation that will allow us to "dec- 
orate" the plaquettes so as to encode the effects of interactions between dimers 
beyond just plaquettes with parallel or non-parallel dimers. To this end, let us begin 
by rewriting the RK Hamiltonian in terms of local objects i^Q^ that denote one of the 



two possible flipppable dimer configurations on plaquette p, i.e., either fff^ = 



or i^Q^ = HJp (more formally, il^^ G {Illp, Dp})- Given the flippable plaquette 

Iq'^ we denote by i'o'^ the flippable plaquette p obtained from a 90° rotation of the 
(p) ^ if p(p) - I-f e(p) _ l-r if p(p) 



p 



dimer covering in (i.e., if t^' = IHp then = Up, or if = Op then 

^ = Up) 

The RK Hamiltonian then becomes 

^iJK = ^E<5,(.) (4) 

where the sum over the local operators 
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runs over all flippable plaquettes. From now on we will drop the explicit depen- 
dence on the plaquette index p in the notation for a flippable plaquette: = £q. 

With this notation in hand, one can extend the definition of a flippable plaquette 
Iq to that of a decorated flippable plaquette 1^ in the following way. Whereas 1^ 
specifies the position of the plaquette and the orientation of the two parallel dimers 
covering it, £q contains additional information on the configuration of dimers at 
neighboring plaquettes. This information can be encoded in a vector m e (0, 1}" 
that lists whether each of n neighboring bonds is occupied by a dimer or not. For ex- 
ample, m = (m^, . . . , 777.4) when the additional information corresponds to speci- 
fying whether the four bonds that face the four edges of the plaquette p are occupied 
by dimers (m^ — 1) or not {m^ — 0), as is illustrated in Fig. 1. Given a flippable 

m, m, 



VOr, 



•P 



mj 



mj 



m3 



Fig. 1. Example of a local update 1^ ^ £q between decorated flippable plaquettes. A 
decorated flippable plaquette is defined by the choice of the plaquette position p on the 
square lattice, by the orientation of the two parallel dimers covering the plaquette p (e.g., 
horizontal for £q and vertical for Iq in this example), and by the values or 1 taken by the 



four parameters m,, i = 1, 



4 defined on the four bonds that face the four edges of 



plaquette p. The value = 1 (mj = 0) corresponds to bond i being covered (not covered) 
by a dimer. In this example m = (1, 0, 0, 0). 

plaquette i*, the operation — > is defined by flipping the flippable plaquette 
p (see Fig. 1). As before, to any decorated flippable plaquette £q corresponds one 
and only one £q. A straightforward generalization of the RK Hamiltonian (4) then 
follows by replacing the sum over flippable plaquettes with the sum over decorated 
flippable plaquettes, 

H^l^Qe*,, (6) 

where 



Q,. := ^(m) \Q{Q + ^;(m) - t(in) + 



(7) 



This Hamiltonian remains local although now the coupling constants v{ni) , v(m) e 
M and t(m) G M can be used to encode interactions that extend beyond the two par- 
allel dimers at a given flippable plaquette and involve, for example, the four dimers 
belonging to the four plaquettes that share a bond with plaquette p in the case of 
Fig. 1. The GS for generic values of v(m), v(m) e R and t(m) e R is not known 
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in closed form. However, when 



v{m)v{m) = r{m), t(m) > 0, (8) 

the GS is given by 

\^) := Y: CiC)\C) (9) 
ce5o 

within any irreducible sector of the Hilbert space under the action of the RK Hamil- 
tonian provided the following integrability condition relating the local data f (m) 
and the global expansion coefficient C{C) holds [33]. For any dimer configuration 
C and for any decorated fiippable plaquette £q present in C the dimer configuration, 
C is uniquely defined by replacing Iq with Iq, while the integrability condition is 
satisfied whenever 

C(C)/CiC) = i (t;(m)/t(m) + t{ni)/^J{m}) (10) 

holds. The conditions above precisely define the Stochastic Matrix Form (SMF) 
decomposition of a quantum Hamiltonian discussed in Sec. 1 . For the remaining of 
this section, we will focus on the example of decorated fiippable plaquettes depicted 
in Fig. 1. 

The existence of the global expansion coefficients in the GS (9) satisfying the inte- 
grability condition (10) can be verified for the choice 



t(m 



m) = exp (u6N^JJ/2T) , 

exp(u6Nil^/2T] , (H) 



with 



SN^P = ±[(mi + mg) - (m2 + m^)] = -5N^\ (12) 

Here, m/T G M and the + (— ) sign is associated to the vertical (horizontal) orienta- 
tion taken by the two parallel dimers occupying the fiippable plaquette p in £q. One 
verifies that SN^P can only assume the values 0, ±1, and ±2. With this choice, the 
Hamiltonian 

is the sum over (non-commuting) operators with an index running over all deco- 
rated fiippable plaquettes 

uSN^V /2T u5N^^/2T_ _ _ _ 

Q,:/-=e \Q{Q+e \to){to\ - \to){il\ - \Q{to\, (14) 

each of which is proportional to a projection operator: 

/ /2T uSn1L''/2T\ _ 

Q,l= e"'^^5/^%e ^5 ' )q,,. (15) 
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The nodeless wavefunction 



^o)= Ee^'^^^IO, (16) 



where A*'^''^'' is the number of flippable plaquettes in the classical dimer configura- 
tion C, is annihilated by the action of each for all as follows from verifying 
that 

6NiP = Ni^^ - N^. (17) 

Therefore, |\1'q) is a GS of Hq, which, with the help of Perron-Frobenius theorem, 
can be shown to be unique within each irreducible sector of H^^^ under the action 



of Hq [|33n . Remarkably, the GS expectation value of any quantum operator O that 
is diagonal in the preferred basis {|C)} can be written in term of an equilibrium 
thermal average for a square lattice classical dimer model. 



(18) 



(^o|0|^o) _ ^ Yr[4'^+N^c'^){C\0\C') 



with 

Z{T/u) := Yl e-^c"VT^ _ _^^(/)^ ^^g-^ 



and 

Oc := (C|0|C). (20) 

It follows that the zero-temperature phase diagram of the interacting quantum SLDM 
(13) contains the phase diagram of the interacting classical SLDM (19). The next 
section is devoted to the numerical study of the phase diagram of the interacting 
classical SLDM (19). 



3 The associated classical model 



The interacting classical SLDM defined by the partition function (19) has been 
extensively studied by Alet et al. in Ref. [ilOi1 for positive values of the coupling 
constant u. Notice that in the range K = u/T ^ (0, oo) the classical energy E^^^ 
favors configurations with a large number of flippable plaquettes, while the diago- 
nal term in the quantum Hamiltonian in Eq. (13,14) always penalizes the presence 
of flippable plaquettes for all values of K\ Conversely, in the range K E (— oo, 0) 
the presence of flippable plaquettes is also penalized at the classical level. Configu- 
rations with the largest number of flippable plaquettes are referred to as the colum- 
nar state (every dimer has two parallel neighboring dimers along every other row 
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or column). A representative among all configurations with no flippable plaquettes 
is the staggered state. A staggered state is obtained from its parent columnar state 
upon translation of every other dimer of each column, say, by one lattice spacing 
along the direction parallel to the dimers. For brevity, we will refer to the parameter 
range K > a.s the columnar side of the interaction, and to the parameter range 
< as the staggered side of the interaction. 

This section is devoted to studying the full parameter range K e (—00,00) of the 
interacting classical SLDM (19) using transfer matrix techniques. A brief summary 
of Alet's results is given in Sec. 3.3.1, where we also present our results on the 
columnar side of the interaction. 



3.1 The conformal field theory description in the T = 00 limit 



At infinite temperature (K = 0), the partition function (19) reduces to the non- 
interacting classical SLDM. It exhibits critical spatial correlation functions with 
power-law decay Jsi 3]. The long-wavelength limit of this model is known to be 
described by the 2D Sine-Gordon field theorywhose action can be written in terms 
of a continuous (height) scalar field h [flsl lslM. 



S= d\ 



7r5(|V/i(r)p + r cos (27rg/i(r 



(21) 

where g is called the stiffness and the V > Q term is called the locking poten- 
tial. Remarkably, the stiffness and the periodicity of the locking potential are fixed 
uniquely to the values g = 1/2 and g = 4, respectively, if the Sine-Gordon ac- 
tion (21) is to encode the long distance asymptotics of the non-interacting SLDM 

. The scaling dimensions „ of the so-called vertex operators at the free-field 
fixed point V = of the Sine-Gordon field theory (21) can be classified in terms of 
their electric and magnetic charge e and m, respectively 



e, m G Z. (22) 



The physical interpretation of electric vertex operators is that their correlation func- 
tions represent the long distance asymptotics of the dimer correlation functions 
provided their charge e is a multiple of g = 4. The physical interpretation of mag- 
netic vertex operators is that their correlation functions represent the long distance 
asymptotics of the monomers correlation functions whereby it is understood that 
monomers are defects in a dimer covering by which sites are not the end points 
of dimers. Monomers will be introduced at the microscopic level in Sec. 5 but are 
absent in the present SLDM, i.e., m = must be enforced. If so, the most relevant 
electric vertex operator with e multiple of g = 4 is the locking potential with scal- 
ing dimension = 16 in the i^T = limit. We conclude that the free-field fixed 
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point V = is the attractive fixed-point of the Sine-Gordon theory (21) if it is to 
capture the long-distance physics of the non-interacting SLDM. 



3.2 Construction of the transfer matrix 



In order to study the phase diagram of the interacting classical SLDM (19) at finite 
values of the reduced coupling constant K — u/T, it is convenient to use a combi- 
nation of numerical transfer-matrix (TM) calculations in the infinite-strip geometry 
with conformal field theory (CFT) arguments. 

We define the interacting classical SLDM (19) on an L x M, M = oo square 
lattice and impose periodic boundary conditions in both directions (i.e., wrapped 
around a torus with infinite principal radius). The width L must be even to respect 
the bipartite nature of the lattice. The TM T^^^ connects one row of the lattice to 
the following one along the principal (infinite) axis of the torus, as illustrated in 
Fig. 2, and satisfies 

Z{T/u; L, M) = Tr [t^^^ (T/u)] ^ , (23) 

where Z{T/u;L, M) is the partition function of the system. We shall assume (for 
simplicity) that the TM T^^^ can be diagonalized through a similarity transforma- 
tion, and we label its (positive) eigenvalues in descending order 

r(^) ^ diag ( a(^) a(^) ) . . . ) , Ai") > Ai^) > a(^) > • • • . (24) 

Anticipating an exponential growth with L of the TM eigenvalues, we also define 
the exponents 

/„(L):=-ilnAW, n = 0,l,2,---. 



L 

The sign is chosen here by convention 



(25) 



The dimensionless intensive free energy in the thermodynamic limit 

/(r/«):= lim/(r/M;L) (26) 

iy— >oo 

is related to the TM T^^) by 

f{T/u-L) -i Jim^^lnTr [r(^)(r/K)]'' (27) 

and is thus solely controlled by the largest (non-degenerate) eigenvalue Aq^^ of the 
TM T^^\ 

f{T/u) = - lim \ In h.f\T/u) = lim f^{T/u- L). (28) 

L—>oo ±j L^oo 
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If the long wavelength limit of the interacting classical SLDM (19) is captured by 
a CFT, the central charge of the CFT and the scaling dimensions of the primary 
fields in the CFT can also be extracted from the finite size (L) dependence of the 
TM T^^\ The dependence on T/u of the central charge c{T/u) is given by 



fiT/u-L) = f{T/u)-^^^^ + Oil/L') 



(29) 



The dependence on T/u of the scaling dimensions d^{T/u) of the CFT primary 
fields is in turn given by 



f^{T/u;L)-f{T/u;L) 



^^^ + 0(1/L3), . = 1,2, 



(30) 



We now turn to the explicit construction of the TM T^^^. To this end, we introduce 
the variables on the bonds of the lattice, = (1) if the edge is empty (oc- 
cupied) by a dimer, as is illustrated in Fig. 2. Any allowed (initial) configuration 



n2L-i 


n'2L-i 
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n2L-4 


■^2L- 2 

n2L-l 





A 





Fig. 2. Labeling of the edge variables in the TM T^^) . Periodic boundary conditions are 
assumed along the x direction, i.e., = ng^,. 



n = {n^, i = 0, . . . , 2L — 1} must then satisfy 

^2a:-l + + ^2x+l < 1> Vx = l, ...,L, 

while the TM T^^^ connects only configurations n and n' that satisfy 



(31) 



(32) 
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We can then write 



(33) 



where the two contributions due to the interaction U^^^{T/u) and V^P^{T/u) (ac- 
counting for horizontal and vertical parallel dimers, respectively), and the contri- 
bution due to the constraint Tj^^]^{0) take the form 



1^ 



x=l 



ui^iT/u) = exp 
v!.^liT/u) = exp 



U ^^2x''^2x+2 + '^2x'^2a; 




2a;+2 



(34) 



rp ^2a;+1^2a;+l 
. x=0 



respectively. 



The computational effort to obtain the eigenvalues of the TM T^^^ can be reduced 
by looking for quantities that are left invariant under the action of the TM. This 
allows to block diagonalize T^^\T/u) and to compute the largest eigenvalues sep- 
arately in each sector. As discussed in Ref. [10] and using the notation defined in 
Fig. 2, one can show that the quantity 



L/2-1 

E 



n 



4j 



n 



4i+2 



L L 



(35) 



is conserved as one proceeds along the y-axis through repeated applications of the 
TM. One can thus use to label the diagonal blocks of T^^^(T/u). Observe that 
any block with ^ corresponds to having \WJ monomers on the same sub- 
lattice at y = — oo, and \ WJ monomers on the opposite sublattice at y = oo, as 
discussed in Ref. [llOll . We will impose the condition that there are no monomers, 
i.e., = 0, when considering the electrical sector of the CFT in isolation. We 
shall assume that, upon ordering the eigenvalues of the TM within each sector ac- 
cording to 



A 







>A\ 



> 



(36) 



Aq^'^'' is the global maximal TM eigenvalue. If so, it is Aq'''^' that controls the free 
energy and central charge of the system following Eq. (28) and Eq. (29) (with 
the possible exception of the transition to the staggered phase, as discussed in 
Sec. 3.3.2). In Sec. 4.1 we show how the introduction of longer-range couplings 
in the dimer model gives rise to new phase transitions that are not entirely captured 
by the W,^. = sector, and one needs to compute the largest eigenvalue from the 
global TM. We shall also define the exponents 



--lnA(^-^) 
111 J i„ , 



n 



0,1,2, 



(37) 
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within each irreducible sector. Equation (30) together with Eq (22) imply that the 
ratio Ao°'^''/A^°'^'' determines the scaling dimension q while the ratio A'q'^^ /A^''^' 
determines the scaling dimension ^. From the values of either of these two scal- 
ing dimensions one obtains the stiffness g of the CFT of Eq. (21), and therefore 
the scaling dimensions of all the operators in the theory. Measuring both of them 
allows for a further check on the reliability of the numerical results, as the product 
of d^ I and di q must remain constant even if g varies (c^q id^Q = 1/4). 



3.3 The transfer matrix results 



Exploiting the sparse nature of the TM due to the constraint-enforcing term T^\{<S), 
we compute its largest eigenvalues using a combination of hashing techniques to 
encode the state space of the system and routines from the ARPACK library [Isill 
based on the implicitly restarted Arnoldi method (or implicitly restarted Lanczos 
method, whenever applicable), which are particularly suitable to handle large sparse 
matrices. Specifically, we compute the largest and first subleading eigenvalues of 
the TM in the = sector, as well as the largest eigenvalue of the TM in the 
= 1 sector. This allows us to obtain, from finite size scaling, the dimensionless 



intensive free energies /q°^ and Jq'^' in the sectors = and = 1, respec- 
tively, the central charge c, and the scaling dimensions d-^ q and dQ ^ of the electric 
and magnetic vertex operators, respectively. We also compute the global largest and 
first subleading eigenvalues of the TM, i.e., independently of the value of W^. We 
can thus verify that the global dimensionless intensive free energy and the central 
charge are indeed controlled by /g^''. 



(■(1) 



For each value of the reduced coupling constant K = u/T and in each block- 
diagonal sector, we consider systems of size L = 6, 8, 10, 12, 14, 16. We then use 
standard finite-size scaling techniques [|52ll to extrapolate the desired quantities. 
Namely, we fitted the computed values 



/(°HL) = -iln(Ar), 



L 
1 

L 
1 

L 



ln(Af'^)) 
ln(A(^'^)) 



-ln(Ar^) 



(38) 



for L = 6, 8, 10, 12, 14, 16, with the scaling forms 

A 



f- 



TTC 



27r(io,i 
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+ jj, 



(39) 
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respectively. The term is introduced to speed up the convergence of the fitting 
routine for the parameters c, di^, and (io,i [l52i1 . Specifically, we first obtain esti- 
mants of the above parameters by fitting data points with Lq < L < i^max' where 
-^max = 16 is the largest strip width that we consider, and Lq is varied from 6 to 
-^max ^ 4 (i^max — 2 in the case of the scaling exponents). The resulting estimants 
-^max)' '^i,o(-^0' -^max)' ^^'^ '^o,i(-^0' -^max) given in Tables 1, 2, and 3. 
We then extrapolated the estimants in the limit Lq — > oo by assuming the power 
law form [i52i1 

^(-^0' -^max) = 2; + kLQ^, (40) 
for each of the three quantities x = c, difi, do^i, and performing an appropriate 
fit of the data in Tables 1, 2, and 3. Whenever the last three estimants in a row 
of Tables 1, 2, and 3 appear in monotonically decreasing order, they are used to 
obtain the constants x, k, and p. When this is not the case, or whenever the power 
p obtained from the fit is too small to produce a reliable extrapolation (namely if 
p < 1, from experience), then the Ising-like value p = 2 is assumed and only the 
last two estimants are considered for the fit to obtain x and k (denoted by a ( )* in 



the Tables 1, 2, and 3) ll52ll . A rough estimate for the error bar can be obtained by 



considering the variation among the estimants [15211 . In the results shown in Figs. 3- 
6, the error bars are computed as the square root of the variance of the last three 
estimants used for the extrapolation. In most cases, the error bars are smaller than 
the symbols used in the plots. 

The results for the central charge and scaling exponents are shown in Figs. 3 and 4, 
for the columnar and staggered side of the interaction respectively. 

The same fitting approach is used to extrapolate the global free energy / from its 
estimants /(Lq, Lj^^^x)' well as f^^^ (in the sector = 0) and f^^^ (in the 
sector = 1). The data are omitted here for brevity, and we show only the results 
for the staggered side in Fig. 6. 



3.3. 1 The columnar side 

On the columnar side of the interaction (Fig. 3), our results are in good agree- 
ment with the ones obtained by Alet et al. . The system remains critical for small 
but finite values of the reduced coupling K = m/T up to some critical value in 
the range (1.5, 1.7) (compare with the critical temperature yMumnar) ^ 0.65 u in 
Ref. [10]) and its long wavelength limit is described by the CFT (21) with continu- 
ously varying stiffness g = 2 do i. In the limit of K 0, one recovers the expected 
value g = 1/2 (not shown) and the cosine term is irrelevant as d^^ = 8/g = 16. 
When the value of K is increased, the stiffness increases monotonically and the 
scaling dimension of the cosine term decreases correspondingly, until this operator 
becomes marginal at T = 7|(coiumnar) ^^^^ _ 2^ g = 4) system undergoes then a 
Kosterlitz-Thouless (KT) transition into the columnar ordered phase, as confirmed 
by the Monte Carlo simulation results presented in Ref. ifioll . 
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Table 1 

List of the estimants c(Lq, ^max) obtained from parabolic fits in for all the values 
of K = u/T considered in this paper. The error bar on the extrapolated value is obtained 
from the variation among the estimants used in the corresponding extrapolation. 



K 


cfe 16) 


cf8 16) 


cflO 16) 


cfl2 16) 


Extrapolation 


1.7 


0.9642 


0.9510 


0.9445 


0.9380 


0.923 ± 0.005* 


1.5 


1.0135 


1.0083 


1.0099 


1.0101 


1.010 ± 0.001 


1.0 


0.9900 


0.9885 


0.9924 


0.9966 


1.006 ± 0.003* 


0.7 


0.9904 


0.9886 


0.9925 


0.9966 


1.006 ± 0.003* 


0.5 


0.9906 


0.9889 


0.9923 


0.9960 


1.004 ± 0.003* 


0.4 


0.9907 


0.9891 


0.9922 


0.9955 


1.003 ± 0.003* 


0.35 


0.9910 


0.9893 


0.9934 


0.9953 


0.998 ± 0.002 


0.0 


0.9947 


0.9919 


0.9951 


0.9979 


1.004 ± 0.002* 


-0.1 


0.9959 


0.9926 


0.9950 


0.9968 


1.001 ± 0.002* 


-0.22 


0.9974 


0.9935 


0.9944 


0.9960 


1.000 ± 0.001* 


-0.35 


0.9985 


0.9938 


0.9952 


0.9985 


1.006 ± 0.002* 


-0.5 


0.9991 


0.9937 


0.9951 


0.9993 


1.009 ± 0.002* 


-0.7 


0.9963 


0.9923 


0.9932 


0.9947 


0.9981 ±0.001* 


-1.0 


0.9507 


0.9716 


0.9808 


0.9864 


1.005 ± 0.006 


-1.12 


0.8998 


0.9401 


0.9617 


0.9741 


1.007 ±0.01 


-1.28 


0.7890 


0.8563 


0.8990 


0.9276 


0.99 ± 0.03* 


-1.5 


0.5719 


0.6472 


0.7047 


0.7491 


0.85 ± 0.04* 


-2.0 


0.2676 


0.2765 


0.2896 


0.2991 


0.321 ± 0.009* 


-3.0 


-0.0118 


-0.0086 


-0.0063 


-0.0068 


-0.008 ± 0.001 


-4.0 


0.0012 


-0.0006 


-0.0001 


-0.0001 


-0.0001 ± 0.0002 



Notice that, contrary to the RK model in Eq. (1) Bol . in our quantum model the 
presence of a possible resonating plaquette phase seems to be excluded according 
to the Monte Carlo results by Alet et al. [|10|1 . 



3.3.2 The staggered side 

On the staggered side of the interaction (Fig. 4), the central charge remains again 
constant at c = 1 for small but finite values of the reduced coupling K, until it 
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Table 2 

List of the estimants di^o{LQ, ^niax) obtained from parabolic fits in for all the values 
of if = u/T considered in this paper. The error bar on the extrapolated value is obtained 
from the variation among the estimants used in the corresponding extrapolation. 



K 


d^ n(8 16) 


di nflO 16") 


di nfl2 16) 


d^ 0(14 16) 


Extrapolation 


1.7 


0.0878 


0.0838 


0.0806 


0.0779 


0.071 ± 0.002* 


1.5 


0.1489 


0.1483 


0.1478 


0.1475 


0.1465 ±0.0003* 


1.0 


0.3095 


0.3090 


0.3087 


0.3084 


0.3076 ± 0.0003* 


0.7 


0.4352 


0.4350 


0.4348 


0.4347 


0.4344 ± 0.0001* 


0.5 


0.5454 


0.5458 


0.5461 


0.5462 


0.5467 ± 0.0002 


0.4 


0.6120 


0.6128 


0.6132 


0.6135 


0.6141 ± 0.0003 


0.35 


0.6488 


0.6498 


0.6504 


0.6507 


0.6517 ±0.0004 


0.0 


0.9986 


0.9993 


0.9995 


0.9997 


0.9998 ± 0.0002 


-0.1 


1.1384 


1.1381 


1.1379 


1.1379 


1.1379 ±0.0001 


-0.22 


1.3375 


1.3362 


1.3359 


1.3358 


1.3358 ± 0.0001 


-0.35 


1.5979 


1.5998 


1.6008 


1.6012 


1.6016 ± 0.0006 


-0.5 


1.9848 


1.9984 


2.0018 


2.0026 


2.003 ± 0.002 


-0.7 


2.1119 


2.1336 


2.1066 


2.0768 


1.99 ±0.02* 


-1 n 








9 41 59 


a n + n 1 
in u. 1 


-1.12 


1.7082 


1.9126 


2.0810 


2.2189 


2.6 ±0.1* 


-1.28 


1.4078 


1.5727 


1.7151 


1.8366 


2.2 ±0.1* 


-1.5 


1.0517 


1.1546 


1.2495 


1.3323 


1.56 ±0.07* 


-2.0 


0.4538 


0.4469 


0.4471 


0.4499 


0.453 ± 0.001 


-3.0 


0.9469 


0.9481 


0.9534 


0.9605 


0.980 ± 0.005* 


-4.0 


1.6293 


1.6312 


1.6326 


1.6340 


1.638 ± 0.001* 



abruptly drops to zero at some critical value in the range (—1.5, —1.28). We can 
therefore estimate the critical temperature for the transition to the staggered phase 
to be T = T^^^^&&^^^) ~ (-0.72 ± 0.05) u. The system remains critical for tempera- 
tures above T^Cstaggered) Yv^hilc the stiffness g decreases monotonically with decreasing 
T, as signaled by the behavior of the scaling dimensions di and (io,i- 



Notice that, on this side of the interaction, all the allowed operators with zero mag- 
netic charge (i.e., those with e = An, n = 1,2,...) have scaling dimensions (ie,o 
larger than two. The absence of an allowed operator that becomes relevant for 
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Table 3 

List of the estimants do,i {Lq , ^niax) obtained from parabolic fits in for all the values 
of K = u/T considered in this paper. The error bar on the extrapolated value is obtained 
from the variation among the estimants used in the corresponding extrapolation. 



K (io,i(8, 16) 4,1(10, 16) do,i(12, 16) do,i(14,16) Extrapolation 



1.7 


1.8139 


1.9017 


1.9671 


2.0186 


2.16 ±0.05* 


1.5 


1.4269 


1.4689 


1.4972 


1.5175 


1.57 ±0.02* 


1.0 


0.8044 


0.8095 


0.8121 


0.8136 


0.817 ±0.002 


0.7 


0.5735 


0.5747 


0.5752 


0.5755 


0.5759 ± 0.0003 


0.5 


0.4562 


0.4567 


0.4569 


0.4570 


0.4572 ± 0.0001 


0.4 


0.4061 


0.4064 


0.4066 


0.4067 


0.4067 ± 0.0001 


0.35 


0.3829 


0.3832 


0.3833 


0.3834 


0.3835 ± 0.0001 


0.0 


0.2494 


0.2497 


0.2498 


0.2498 


0.2499 ± 0.0001 


-0.1 


0.2191 


0.2194 


0.2195 


0.2196 


0.2196 ± 0.0001 


-0.22 


0.1865 


0.1868 


0.1869 


0.1870 


0.1871 ±0.0001 


-0.35 


0.1554 


0.1557 


0.1558 


0.1559 


0.1559 ±0.0001 


-0.5 


0.1243 


0.1246 


0.1248 


0.1248 


0.1249 ±0.0001 


-0.7 


0.0901 


0.0902 


0.0904 


0.0905 


0.0908 ± 0.0001 


-1.0 


0.05186 


0.05166 


0.05164 


0.05164 


0.05164 ±0.00001 


-1.12 


0.0404 


0.0400 


0.0399 


0.0399 


0.0398 ± 0.0001 


-1.28 


0.0279 


0.0275 


0.0273 


0.0271 


0.0266 ± 0.0002* 


-1.5 


0.0089 


0.0110 


0.0147 


0.0145 


0.014 ± 0.002* 


-2.0 


-0.0099 


-0.0005 


-0.0052 


-0.0042 


-0.001 ± 0.002 


-3.0 


-0.0011 


-0.0003 


0.0002 


0.0002 


0.0001 ± 0.0003 


-4.0 


0.00008 


0.00002 


0.00000 


0.00000 


0.00000 ± 0.00001 



T > j'j'^'aggered) conslstcnt with the transition being first order, as expected in 
analogy with the RK Hamiltonian 1I12I1 (see also the results on the free energy pre- 
sented below). 



Finally, the stiffness g seems to vanish exactly at the staggered transition (within 
numerical error), suggesting that a tilting transition occurs at j'^(staggered) briefly 
noted in Ref. (w]- Observe, however, that our estimate T^^'^ss^'^'^) ~ (-0.72 ± 
0.05)m differs from the estimate Tj'*'^gg'=^'='^' ~ -0.449(1)m in Ref. |fl^. A possible 
explanation for this discrepancy might be due to strong topology-dependent finite- 
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1/K = T/u 



Fig. 3. Central charge and scaling exponents for the electric and magnetic monopole op- 
erators on the columnar side of the interaction. Following the convention in Ref. [10], we 
plot the scaling exponents = 2diQ and = 2 ^ (equal to l/g and g from Eq. (22), 
respectively), instead of the scaling dimensions ^. The results by Alet et al. are also 
shown for comparison iiod . 




1/K = T/u 



Fig. 4. Central charge and scaling exponents for the electric and magnetic monopole oper- 
ators on the staggered side of the interaction. Following the convention used in Ref. [101 . 
we plot the scaling exponents = 2d^Q and = 2 ^ (equal to l/g and g from 
Eq. (22), respectively), instead of the scaling dimensions ^. The divergence of close 
to T/u ~ —1 makes it increasingly difficult to obtain a reliable extrapolation for L — > oo, 
hence the large fluctuations observed. The smaller parameter Um is much less affected and 
the curve 1/am can be used as a guide to the eye for the diverging behavior of a^. 

size effects. For example, the central charge of the system computed using the 
global transfer matrix as is done in Fig. 5 suggests a drop of the central charge 
at the much lower value of —0.3m than the value —0.72m extracted from Fig. 4. In 
turn, this suggests a strong dependency on the topological sectors of the finite-size 
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estimate for the critical temperature. Understanding this interplay between finite- 
size effects and topological sectors is beyond the scope of this paper. 
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Fig. 5. Estimants of the central charge and scaling exponent for the magnetic monopole 
operator on the staggered side of the interaction, derived using the global transfer matrix 
for L = 8, 12, 16 [53]. Following the convention used in Ref. I IQ] , we plot the scaling 



exponents 
scaling dimensions 



2 ^iid '^r, 



2dQi (equal to 1/g' and g from Eq. (22), instead of the 



For comparison, we also show the estimants of 



2d 



corresponding to the most relevant operator in the CFT, obtained from the first subleading 
eigenvalue of the global TM. As expected, a„ = a^j^ as long as the CG description in 
Sec. 3.1 holds, i.e., for \T/u\ > \{T / uf^^^''''''^. 



3.3.3 The free energy 

The results for the global free energy of the system (/) and for the constrained free 
energies in the = (/o°^) and = 1 (fo^^) sectors are shown in Fig. 6 for 
the staggered side of the interaction. The continuous KT transition on the columnar 
side is in general difficult to detect from free energy measurements, and the curves 
(not shown) are indeed smooth and featureless across t^^°^'^^'^^^'> _ 

On both sides of the interaction and for the range of temperatures considered here, 
the three free energies are indistinguishable (/ = flf^ = /q^'*). This overlap is 
indeed expected since the number of sectors is linear in system size, and the 
presence of defects at y = ±oo should not affect significantly the entropy of each 
sector. 

As mentioned above, the transition on the staggered side is expected to be first 
order. In fact, not only has the set of all possible staggered configurations vanish- 
ing entropy, but they also do not allow for local fluctuations (without violating the 
dimer constraint), so that even the set of all the configurations connected to a stag- 
gered one through a single thermal fluctuation has vanishing entropy. This peculiar 
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-10 -8 -6 -4 -2 

1/K = T/u 



Fig. 6. Free energy of the system obtained from the global TM and from the TM restricted 
to the = and = 1 sectors, respectively, for the staggered side of the interaction. 
The inset shows a close-up view of the temperature range where the transition is expected to 
happen, and the dashed line indicates the value of j'^^^^s^^^^'^^ predicted from CFT arguments 
in Sec. 3.3.2. 



connectivity of phase space close to the staggered phase leaves little room for a 
possible continuous phase transition, and indeed in the RK Hamiltonian (1) the 
transition to the staggered phase is similarly observed to be first order. 



The accuracy in our data does not allow for a precise identification of the transition 
temperature from the free energy plot. However, the overall behavior is in agree- 
ment with the value of T^^^itaggered) obtained from CFT considerations, as shown in 
the inset in Fig. 6 (notice that the staggered phase has both vanishing entropy and 
vanishing energy, therefore a lower bound for the transition temperature is obtained 
from the location of the point where the free energy first vanishes). A deeper in- 
sight on the nature of this transition could perhaps be obtained from Monte Carlo 
simulations, even though the non-local nature of the thermal fluctuations above the 
staggered phase are likely to lead to dynamical slowing down and glassiness, as 
observed in similar coloring models [|54|] . A variational (cluster) mean field ap- 
proach nSSL 15411 may prove more powerful in obtaining an accurate estimate of the 
transition temperature to compare with the one from CFT arguments. Such analysis 
is however beyond the scope of the present paper. 



An intriguing question about this transition to the staggered phase would be to 
understand what happens in this model to the so-called devil's staircase scenario, 
proposed in the original Rokhsar-Kivelson square lattice dimer model (1) [|26|l. 
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4 Stability of the critical line 



The question that we want to address in this section is how stable is the critical 
segment —1/0.72 < u/T < 1/0.65 in the zero-temperature phase diagram of the 
interacting quantum SLDM (13) to a perturbation by some longer-range interaction 
E^^'' between dimers. In order to address this question, at least from a qualitative 
point of view, we need to substitute 

4") ^ + (41) 

in the classical partition function (19) by retracing all the steps that lead from the 
interacting quantum SLDM (13)_to the classical partition function (19). (That this 



is possible was proven in Ref. 113311 .) For the sake of concreteness, let us consider 
the case where the new contribution to the classical energy of a dimer configuration 
C is given by 



= E E -1? hl^ + -S? + -S? + + -S? + -S] ■ (42) 



a=h,v 



Here, the sum over i runs over all lattice sites, the sum over a = h,v runs over the 
occupation number by horizontal and vertical dimers respectively, and n\l^ = 1 (0), 
r = 0, ■ ■ ■ , 6 if the bonds depicted in Fig. 7 are covered (not covered) by a dimer. 
Notice that the J interaction is comprised of a nearest-neighbor part 

(0) (2) , (0) (5) 

>i ^ ' >i ^ ' _l_ (in ^ ' (in ^ ' 

"ia "ia "ia "ia 

and of a next-nearest-neighbor part 

/in (in ^ ^ I (in ,in I (j-n /in I /in (in 

"ia "ia "ia "ia "ia "ia "ia "ia 

that are equally weighted. Therefore, varying both u - which corresponds to a pure 
nearest-neighbor interaction between dimers - and J is equivalent to tuning sep- 
arately the nearest-neighbor and next-nearest-neighbor couplings between dimers. 
In the limit J/T ^ \u\/T, the coupling J favors dimer configurations where the 
dimers are mostly aligned along the same direction. In this limit, the disordered 
critical phase between the two ordered phases is penalized. For this reason we con- 
jecture that the effect of an increasing value of J/T > is to shrink continuously 
the size of the critical segment — 1/0. 72 < u/T < 1/0.65 at J = until the critical 
value (J/T)^ is reached at which the critical line terminates into a tricritical point. 
A similar behavior is observed in the classical phase diagram of the three-coloring 
model in presence of an Ising interaction and of a uniform magnetic field jssll . This 
scenario is depicted in Fig. 8. Beyond this point, i.e., for J/T > (J/T)^, we con- 
jecture that the columnar and staggered phases are separated by a first-order phase 
transition. The slope of the line of first-order transition points for J/T > ( J/T) c is 
positive because the J > coupling favors the staggered phase over the columnar 
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Fig. 7. Given any site i of the square lattice, the locations of the horizontal bond 



variables n 



n 



n 



i,h 
(0) 



(0) (6) 



and the locations of the vertical bond variables 



(6) 

' "'l,V 



are defined in the left and right panel, respectively. Observe that 
the right panel is obtained from the left panel after a 7r/2 counterclockwise rotation fol- 
lowed by the substitution h ^ v. Given the lattice site z, r E {0, • ■ ■ , 6}, and a E {h, v} 

(r) 

the bond variable n]: ^ takes the value 1 if a horizontal (q = h) or vertical (a = v) dimer 

(r) 

is present while n] ^ = otherwise. The contribution from site i to the classical energy 

E^'^^ defined in Eq. (42) is — J for the dimer covering C in the left panel and —3 J/2 for the 
dimer covering C in the right panel. 
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Fig. 8. Phase diagram of the dimer model with both plaquette counting (u) and short-range 
aligning (J) interactions. The dashed line delimiting the staggered phase is expected to be 
a coexistence line (a first-order phase transition), while the solid line between the columnar 
and the disordered critical phase is expected to be a line of KT (infinite order) transition 
points. The quantitative details shown here are derived from numerical and analytical argu- 
ments in Sec. 4.1. Notice that the phase diagram of a classical system with fixed values of 
u and J as a function of temperature is described by a straight line that passes by the origin 
of the axes. 

phase. The point P^, where the continuous KT phase transition to the columnar 
phase meets the first order phase transition to the staggered phase, is reminiscent 
of the RK critical point (2) in the RK Hamiltonian (1) in that the transition is first 
order on one side and continuous on the other side [ISqI . 



These conjectures have been tested numerically with the help of transfer matrix 
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calculations of the central charge and scaling exponents for the system in presence 
of the additional energy term Ej.^^ in Eq. (42). The numerical techniques that we 
shall use follow step by step those described in Sec. 3, provided we introduce two 
extra factors U^^i{T/J) and vi^liT/J) in Eq. (33), 



U^hT/J) = exp 



x=0 ^ 



((1 - ^2x-\ + %x-2 + ^2x-Z) + ^2x-2 + ^2x-2 
+ (1 — n2x+l + ^2x+2 + ^2x4-3) + '^2a;+2 + '^2a;+2^ 



y^^l{TlJ) = exp 



rp ^2x- + l \ "'2x-l 

. a:=0 



riL-i + n' 



2x+l 



+ n'. 



2x+3 



(43) 



for horizontal and vertical dimers respectively. The notation used here for the la- 
beling of the bonds in the lattice is the same as in Sec. 3.2 (see Fig. 2). The results 
obtained from these numerical calculations are summarized in Sec. 4. 1 and are in 
agreement with the conjectured phase diagram in Fig. 8. 



4.1 Numerical results 



We ran transfer matrix simulations at fixed values of J/T = —0.2, 0.2, 0.5, 1.0, 
and 1.5, while varying u/T along the same set of values used for the J = case. 
We also sampled the line u— J{= 1) in parameter space, at different (inverse) tem- 
peratures = 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.7, 1.0, 1.2, 1.5, 1.7, 2.0, 2.5, 3.0. 

Before presenting the numerical results, let us discuss the low-temperature asymp- 
totics of the system for which analytical arguments can be formulated. For the sake 
of simplicity we will discuss here the case of an L x L square lattice, which one 
can then generalize to the M x L case appropriate for TM calculations. As before, 
L is chosen even when imposing periodic boundary conditions so as to maintain 
bipartiteness of the square lattice. 

We begin by considering the case of u/T <^ —1 and J/T = when the system 
lies deep within the staggered phase. The degeneracy of the classical ground state 
manifold in the staggered phase is subextensive, in the sense that the number of 
degenerate configurations is exponential in the linear size of the system. To see 
this start from the staggered configuration defined in the left panel of Fig. 9(a). 
Draw a line that crosses the mid-points of the occupied bonds along a diagonal, 
as illustrated in the figure. Given the staggered configuration in the left panel of 
Fig. 9(a), and imposing periodic boundary conditions, there are L/2 such diagonals 
(recall that L must be even to ensure bipartiteness). For any such diagonal line, 
we are free to take the dimers that intersect the diagonal and rotate them by 7r/2 
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(a) Staggered phase degeneracy for J = 
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(b) Columnar phase degeneracy for u = J 

Fig. 9. Examples of ordered phases in the asymptotic limits n/T— >— cxd,J = (Top) and 
u/T,J/T — > oo (Bottom). For J = all staggered configurations are degenerate, and can 
be obtained from the four canonical configurations - one of which is shown in the Top-left 
figure - by 7r/2-degree rotation of the dimers along diagonals (Top-right). Equivalently, for 
u = J all configurations obtained from one of the four columnar states by one-sublattice 
translations of rows of coUinear dimers (Bottom-left ^ Bottom-right) are also degenerate. 

without violating the dimer constraint or creating flippable plaquettes, as it has 
been done on the right panel of Fig. 9(a). Hence, both configurations depicted in 
Fig. 9(a) are degenerate in energy when J /T = 0. By global n/2 rotations of all 
the bonds intersecting any one of the L/2 diagonals, one generates a total of 2^/^ 
degenerate configurations when J/T = [57]. This degeneracy equals the number 
of all possible Ising spin configurations along a chain made of L/2 sites. Notice that 
we could have equivalently well chosen a top-left to bottom-right diagonal instead 
of a top-right to bottom-left diagonal as illustrated in Fig. 9(a). This choice can be 
shown to generate another 2^/^ degenerate configurations, only one of which has 
already been counted before (for a total of 2 x 2^/^ — 1). Similarly, one can show 
that if we take as the starting configuration the one that obtains from the left panel 
of Fig. 9(a) by means of a translation by one lattice spacing in the horizontal (or 
vertical) direction, the same procedure would yield another 2 x 2^/^ — 3 independent 
degenerate configurations (for a grand total of 4 x 2^/^ — 4). 

Switching on the interaction J in Eq. (42) lifts this macroscopic degeneracy by 
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assigning an energy 



=4"'"^ + 4"'^^ (44) 

L-stagg ^stagg *-'Stagg ^ ^ 



to each configuration in this manifold. The energy e'q '^ has been decomposed into 
two contributions. The first one, 

- X (L/2) 

o i/2 (45) 

does not lift the degeneracy of this manifold and can be thought of as resulting 
from the zero-point energy of a diagonal labeled by the integer % — 1, • • • , L/2. 
The second one, 




can be thought of as resulting from an Ising interaction between consecutive diag- 
onals whereby the Ising spin 8,^ variable takes the value +1 (—1) if the diagonal 
intersects vertical (horizontal) dimers. We have thus mapped the problem of eval- 
uating the energy (42) for all configurations obtained from the seed configuration 
in Fig. 9(a) by flipping rigidly all the dimers that intersect top-right to bottom-left 
diagonals as in Fig. 9(a) into an effective one-dimensional Ising model with the 
extensive exchange energy J L/2. Observe that the fact that the nearest-neighbor 
Ising interaction is proportional to the system size in the ID Ising model (46) im- 
plies that (a) the system orders for any non-vanishing value of the dimensionless 
ratio J/T and (b) excitations are separated from the ordered state by an energy gap 
that grows linearly with L. The Ising order is ferromagnetic when J > and an- 
tiferromagnetic when J < 0. The dimers are all aligned along the same direction 
when J/T > (parallel staggered order). When J/T < 0, the dimers are aligned 
along the same direction within a diagonal but are rotated by 7r/2 between any two 
consecutive diagonals (alternating staggered order). Equivalent results are obtained 
if we consider top-left to bottom-right diagonals or a translated starting onfigura- 
tion, as discussed above. Therefore, the ground state manifold of energy (42) is 
four-fold degenerate for any finite J. 

Next, we consider the case of u/T ^ 1 and J/T — when the system lies deep 
in the columnar phase. There are four degenerate columnar ground states one of 
which is depicted in the left panel of Fig. 9(b) and the other three are obtained upon 
rotating globally all the dimers by 7r/2 or upon a global translation by one lattice 
spacing. This degeneracy is preserved when switching on the interaction J, which 
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results in the energy 



411 = -^>4i, (forJ>0). (47) 

The equality holds only for the four staggered configurations with alternating dimer 
directions along any two consecutive diagonals. 

Third, we want to estimate the average energy E^^^] of configurations in the disor- 
dered phase due to the coupling J. This is given by 

4? = (2 X 2 X C„„ + 2 X 4 X C„„J L^ (48) 

where C„„ is the probability to find a pair of parallel nearest-neighbor dimers and 
is the probability to find a pair of parallel next-nearest-neighbor dimers in 
the disordered phase. In the non-interacting limit u, J 0, the probabilities C^^ 
and C„j,„ can be obtained exactly in the thermodynamic limit, as illustrated by 
M. E. Fisher and J. Stephenson in Ref. ^ 

C„n = g, C„nn =4(2"^)- ^^^^ 



Consequently, to leading order in L, one obtains 



2L2 4L2 /I 1 
—r- + 



4 V2 TT 



JL^ (37r - 4) 



(50) 



in the limit u,J ^0. Notice that if J > then ^i^^^ > E^^ > E^/^ . For finite val- 
ues of the couplings J and u, the probabilities C„„ and C„^^ in the disordered phase 
are modified and the above result no longer holds exactly. However, we conjecture 
that the inequality E^fJ > E^^2i,^ holds true throughout the line of critical points 
separating the columnar phase from the staggered phase for J > 0. This conjecture 
is indeed confirmed by our numerical transfer matrix results (e.g., see Fig. 10) that 
clearly indicate a shrinking of the disordered phase for positive values of J. 

For sufficiently large values of J/T, the disordered phase is expected to disappear 
as J favors the parallel alignment of nearest-neighbor and next-nearest-neighbor 
dimers. The transition between the staggered and columnar phases is then con- 
trolled by the competition between the dimensionless free energies of the two or- 
dered phases 

^coi(T) = -^-^-ln(4), 

2^ (51) 
i^stagg(T) = -^-ln(4), 

where for the staggered case we used the energy and entropy consistent with the 
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J > case. This yields an estimate for the transition to occur when J = u that is 
in remarkably good agreement with the numerical results. 

Let us now turn to the TM calculations. When J > 0, the = sector captures 
the physics of the global transfer matrix much like in the J = case and we can 
thus use the same type of approach. The critical nature of the disordered phase ap- 
pears to be preserved by the J coupling, although its width is gradually reduced in 
agreement with our conjectured phase diagram in Fig. 8. For brevity, we show in 
Fig. 10 only the central charge and scaling exponents obtained for J/T = 0.5. The 




K=u/T 



Fig. 10. Central charge and scaling exponents for the electric and magnetic monopole 
operators as a function of K = u/T. Positive values of K correspond to the columnar 
side of the interaction, while negative ones correspond to the staggered side. Following the 
convention in Ref. flO'], we plot the scaling exponents a^ = 2d^Q and = 2 ^ (equal 
to 1/(7 and g from Eq. (22), respectively), instead of the scaling dimensions ^. The value 
of the reduced coupling constant J/T = 0.5 is held constant throughout. The dashed lines 
indicate the location of the staggered and columnar transitions in the J = case discussed 
in Sec. 3.3. 

origin of the peak in the central charge appearing at the first order transition to the 
staggered phase is as of now unclear, though it is likely due to some spurious finite 
size effect (see also the discussion below on the first order transition between the 
columnar and staggered phase at low temperatures). For values of the reduced cou- 
pling J/T larger than (J/T)^ ~ 0.6 the critical region shrinks to the point J = u 
at which a direct transition between the staggered and the columnar phase takes 
place. This is in agreement with the free energy arguments of Eq. (51), suggesting 
that the transition arising from the collapse of the line of critical points is indeed a 
first order transition. Observe that the first-order critical point J = u that appears 
at low temperatures (J/T > 0.6) is special from an entropic point of view. Indeed, 
all the configurations obtained from a columnar one through any combination of 
the translations illustrated in Fig. 9(b) are degenerate in energy when J = u. While 
the number of degenerate lowest-energy configurations is finite in the parallel stag- 
gered and columnar phase, it becomes exponential in the linear size of the system 
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at the J = u critical point. Namely, there are 2^ + 2^ degenerate configurations 
in an L X M system, whereby 2^ (2^) correspond to having all the dimers aligned 
horizontally (vertically). Consequently, at the coexistence point J = m the intensive 
entropy of the system scales with system size as 



In ( 



LM 



(52) 



In turn, this gives a contribution to the dimensionless intensive free energy obtained 
from TM calculations of the form 



In e^ + e*M i 

lim ^---r ^ = -. (53) 

M^oo LM L 



Equation (53) introduces a 1/L finite-size scaling correction, the signature of which 
is to give rise to a peak in the central charge dX J = u{c would actually be infinite 
if the L — > cxD extrapolation could be done exactly) when using the finite size 
scaling Ansatz for the central charge as in Sec. 3.3. We can thus take advantage of 
this effect, when fitting the free energy as in Eqs. (39) and (40), to interpret the peak 
in the free energy when J/T > 0.6 as a signature of a first order phase transition. 
A remnant of the 1/L correction to scaling might also be responsible for the central 
charge peak at the staggered transition that appears in the numerical results when 
J/T < 0.6 (see Fig. 10). 

When J < 0, the numerical analysis becomes more delicate because the stag- 
gered order that the system develops at large, negative values oiu/T has alternat- 
ing dimer directions between consecutive diagonals [see Fig 9(a)]. Studying this 
phase demands larger lattices due to the symmetry breaking pattern which allows 
only values of L that are multiple of 4 in the notation used in Fig. 2. Moreover, this 
type of order lies in the — ±L/4 sector and the transition cannot be captured 
by considering the = sector alone (which is the correct one for the columnar 
transition). In Fig. 11 we present the results obtained from the L = 8, 12, 16 data 
with the global TM. The small number of accessible system sizes does not allow 
for an appropriate L — > cxd extrapolation (40) and the values in Fig. 1 1 are directly 
obtained from the finite size scaling fits (39). As expected, upon changing the sign 
of the coupling J the disordered phase is now favored over the two ordered ones 
and the line of critical points appears to expand with respect to the J = case in 
Fig. 1 1 (Top). From the behavior of the global and restricted free energies one sees 
that the correct asymptotic behavior (47) lim„/r^_oo(/) = — J/2 = 0.1 is captured 
only by the global TM. 
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Fig. 11. (Top) Central charge and scaling exponents for the electric and magnetic monopole 
operators as a function of K = u/T. Positive values of K correspond to the columnar 
side of the interaction, while negative ones correspond to the staggered side. Following 
the convention in Ref. I.lft1 . we plot the scaling exponents = ^d^^ and a„ = 2dQi 
(equal to 1/g and g from Eq. (22), respectively), instead of the scaling dimensions dgm- 
The value of the reduced coupling constant J/T = —0.2 is held constant throughout. The 
dashed lines indicate the location of the staggered and columnar transitions in the J = 
case discussed in Sec. 3.3. (Bottom) Free energy of the system obtained from the global 
TM and from the TM restricted to the W^^ = and W^^ = 1 sectors respectively. Notice 
that only the global free energy saturates at — J/2 = 0.1 for large negative values of u/T, 
as expected from Eq.( 47). 



5 The role of defects 



The purpose of this section is to generalize Hamiltonian (13) in such a way that the 
dynamics of defects becomes possible without spoiling the existence of a GS of the 
form (16). We shall only be concerned with point defects that we call monomers. 
A monomer is a site of the square lattice that is not the end point of a dimer. 
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Fig. 12. There are eight ways for a plaquette of the square lattice to accommodate a single 
monomer (a site occupied by a filled circle) and a single dimer on any of its two remaining 
free edges. These eight configurations are grouped into pairs that define the local moves in 
configuration space. Any of these local moves amounts to a rotation by 7r/2 of the dimer 
and a diagonal hop of the monomer. 

m^ mj 



m^ 



mj 



mj 



Fig. 13. One of the two pairs of decorated plaquettes £1 and £1 that involve hopping of the 
monomer between the upper left and lower right comers of the plaquette p. 

The strategy that we will follow consists of three steps. First, we define the enlarged 
classical configuration space S describing a dilute SLDM where each site belongs 
to at most one dimer. Second, we define the classical configuration energy 



(54) 



(f) 

where Nq ' is, as before, the total number of plaquettes occupied by two parallel 
dimers in configuration C while is the total number of monomers in configu- 
ration C. The energy scale /i G M thus plays the role of a chemical potential for 
the monomers. Third, we need to generalize the local moves fg <^ in the classi- 
cal configuration space Sq that we encountered in Fig. 1 to endow monomers with 
quantum dynamics. To this end, we first define decorated plaquettes that involve 
monomers. Any plaquette of the square lattice has four ways to accommodate a 
single monomer. Given a plaquette occupied by a single monomer, there are two 
ways to accommodate a dimer on one of its two remaining free edges. Any plaque- 
tte of the square lattice has thus eight ways to accommodate a single monomer and 
a single dimer. In turn, these eight ways can be grouped into pairs each of which 
defines a local move in the classical configuration space as shown in Fig. 12. A 
decorated plaquette t\ at p built around a single monomer and a single dimer is 
obtained by taking Fig. 1 and replacing the flippable plaquette p by any of the pla- 
quettes in Fig. 12 as is done in Fig. 13. Then, we want to account for the creation 
and aimihilation of pairs of monomers. Following the same steps as above, we in- 
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(h) 



m, 



m 



(h) 



m 



(V) 

mi 



tl] ^{ fl) 

Fig. 14. A pair of decorated horizontal bonds ^2 and£2 and a pair of decorated vertical 

bonds ^2^''^ and ^^''^ that are related by the process of creation / annihilation of a pair of 
monomers. The information encoded by the decorated bond is the location of the bond, 
whether it is occupied by a dimer or two monomers, and whether the two bonds facing it 
are occupied by a dimer {rn(^2^ = 1) or not (m^^2^ = 0). 

troduce the decorated horizontal and vertical bonds £l^^^ and £l^^\ as well as the 

corresponding local updates £l^^^ ^ £l^^^ and £l^^^ <-> £l^^^ defined in Fig. 14. For 
notational simplicity we will omit in the sequel the reference to the bond orienta- 
tion /i, V, as we have already done for the bond b and plaquette p position labels, 
and for the eight different types of decorated monomer-dimer plaquettes. 

The new quantum Hamiltonian can then be assembled from three separate parts, 
the pure dimer contribution Hq and two terms arising from the introduction of 
monomers. Hi and H2, 



(V) 



m^ 



(V) 



m, 



(V) 



m^ 



-f^tot + Hi + H2 



4? 
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(55) 

( f) 

The energy scales ^ ^1 ^ ^iid t2 are positive. For any n G {0 , 1 } , the integer 6N^: = 

( f) ( f) ( f) 

—6Nj^ is nothing but the change N)f — N^, ' in the number of flippable plaquettes 

induced by the local update ^ when the action of |^*)(^* | on the state \C) 
is non- vanishing (or, equivalently, when the action of (£* | on the updated state 
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\C) is non-vanishing). The integer ^Af^, = —5M-^ takes the value +2 when two 
monomers are created and the value —2 when two monomers are annihilated under 
the local update £3 ^ ^2- We leave it as an exercise to the reader to express 5N\P 

in terms of the integers 2 3 4 defined in Fig. 13 and in terms of the integers mi 2^ 
defined in Fig. 14. 

As desired, the GS of Hamiltonian (55) is given by 

l^t^t) = ^ e(«iv(«+M/c)/2T|c) (56) 

CG5 

where the summation has been extended to account for the enlarged Hilbert space. 
By construction monomers always occur in pairs with one half of the monomers 
residing on one of the sublattice of the square lattice. The zero-temperature phase 
diagram of the quantum Hamiltonian (55) thus contains the finite-T phase diagram 
of the classical system with the partition function 

Z{T/u,fi/T) := 2^ exp I — I = exp I I . (57) 



A detailed study of the phase diagram of this partition function is beyond the scope 
of the present paper. However, according to some preliminary results by Alet et al. 



in Ref. [110], we anticipate a rich and interesting structure, qualitatively illustrated 
in Fig. 15. Notice that the u = 0, yU = —00 limit realizes precisely the RK point (2). 
Also, although the effect of a finite concentration of static monomers on the clas- 
sical critical phase is an interesting open problem, one can reasonably expect that 
the spatial dimer correlations will be preserved below a cutoff length scale dictated 
by the average monomer-monomer separation. Therefore, the limit = real- 
izes a classical neutral gas of charged hard-core particles (the monomers) at a fixed 
density, coupled to a quasi-critical environment (the dimers). 

Finally, having constructed a quantum dilute SLDM Hamiltonian using the Stochas- 
tic Matrix Form decomposition (55) offers the advantage of obtaining (static) GS 
correlation functions directly from the associated classical system, as pointed out 
in Sec. 1. The possibility to use classical numerical techniques, such as Monte 
Carlo simulations and transfer matrix calculations, in the same number of dimen- 
sions gives access to much larger system sizes than quantum techniques, such as 
quantum Monte Carlo or exact diagonalization routines, do. These classical results 
allow to contrast quantum dynamical correlation functions of monomers against 
static correlation functions inherited from the classical partition function (57). In 
this context, related quantum dimer models with mobile holes have been studied 
by exact-diagonalization. Green-function, and classical Monte Carlo techniques by 



Poilblanc and collaborators 115 811 . 
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tricritical point 



dimer + monomer 
liquid 



1 

0.58 



staggered phase 




weakly doped 
columnar phase 



X u/T 

columnar phase 



RK point 

line of critical points kT transition 

Fig. 15. Conjectured phase diagram of the classical partition function (57) after Alet et 
al. in Ref. jlOll that is inherited by the SMF quantum dimer model (56). While for 
u/T < ^^/^"^(coiumnar) jj^g introduction of monomers is a relevant perturbation - which de- 
stroys the critical phase - it becomes marginal at jXcoiumnar) Therefore, a phase boundary 
//(T) is expected to depart from the KT transition point at (T, //) = (r^(coiumiiar) ^ -oo). In 
the close vicinity of (T, ^) = (^"^(coiumnar) ^ — oo), this phase boundary separates a weakly 
doped columnar ordered phase from a dilute (non-critical) dimer liquid and the transition 
between these two phases is expected to be continuous (solid line at the phase boundary). 
For both — T//i and u/T large, the phase boundary /i(r) separates phases of matter with 
large gaps and the transition between these two phases is expected to become first order 
(dashed line at the phase boundary). A tricritical point at some finite {T^,fi^) must then 
separate the first order from the continuous behavior along the phase boundary. Alet et al. 
give an estimate for the tricritical point at T^/u = 0.39(4) and fi^/u = 0.25(10) 1.1(^1 . 



5. 1 The non-interacting monomer gas in a critical background 



In closing, we would like to discuss briefly a special case of the Hamiltonian in 
Eq. (55), namely the one we obtain upon choosing = t^^ = and m = 0. Given 
this choice, ^2 becomes just a scale factor and can be conveniently set to 1. Also, 
the information about the integers defined in Fig. 14 is no longer needed and 

can be dropped. The new Hamiltonian thus obtained can be written in a symbolic 
notation along the line of Eq. (1) as 



(58) 

where the summation is over all vertical and horizontal bonds b and the Hilbert 
space is the same as for Eq. (55). The exact GS of this Hamiltonian 

l^.on) = E e^^y'^'lC) (59) 

CG5 
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is then captured by the classical dilute dimer model in the sole presence of a chem- 
ical potential /i for the monomers, 

Z^o. = E e^''^^^- (60) 

In the limit of fx/T — oo, the monomers behave like deconfined virtual parti- 
cles in a critically correlated dimer background, and the GS wavefunction becomes 
exponentially close to the GS wavefunction of the RK Hamiltonian (1) at the so 
called RK point (3). 



6 Conclusions 

In this paper, we constructed an interacting quantum dimer model with a Hamilto- 
nian obeying a Stochastic Matrix Form decomposition f33ll . This allowed us to de- 
fine an associated classical system whose thermodynamic equilibrium correlation 
functions reflect the behavior of some equal-time ground state expectation values. 
Using this correspondence, we showed that the quantum system exhibits a line of 
quantum critical points that separates two ordered GS of the valence bond crystal 
type. One end point of this line of critical points corresponds to KT transition in 
the associated classical system. The other one is a first order phase transition. The 
locality of our SMF quantum Hamiltonian can be used to argue that the dynami- 
cal exponent z is finite, and likely to be 2; = 2, along the line of critical points. If 
so, one can deduce the power-law decay of quantum correlations in time from the 
behavior of the equal-time spatial correlations. In this sense, the quantum system 
inherits a KT transition from its associated classical system. 

We also studied the robustness of this line of critical points to the introduction of 
longer-range - yet finite - dimer-dimer interactions. In particular we showed how, 
for some choice of the additional interaction, one can tune the width of the line of 
critical points until it gradually collapses onto a single, tricritical point, where the 
KT and first-order transition meet. 

Eventually, we included monomers in the quantum system, i.e., we allowed for 
sites that are not occupied by a dimer. In doing so, we preserved the SMF struc- 
ture of the Hamiltonian, thus being able to study the GS properties of the quan- 
tum system via an associated classical dilute dimer model. Using this correspon- 
dence, we derived qualitatively the zero-temperature phase diagram of the quantum 
dilute dimer model, opening the possibility to obtain equal-time monomer corre- 
lation functions from classical numerical techniques. While this allows access to 
much larger system sizes to study zero-temperature spatial correlation functions, 
unequal-time monomer correlations remain the prerogative of quantum MC or ex- 
act diagonalization studies. 
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